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Abstract. Multi-epoch GHz Very Long Baseline Interferometry (VLBI) data of the radio-loud quasar JllOl+7225 were ana- 
lyzed to estimate the proper motion of extended optically thin jet components. Two components separated from the core could 
be mapped at 1.66 GHz, which is consistent with earlier observations. In one case we found evidence of high apparent super- 
luminal motion (fiapp = 22.5 ± 4) at large (deprojected) distances to the core (22 mas ~ 4kpc, atz = 1.46). Typically in other 
quasars such high separation velocities are only found much closer to the core component. Furthermore the Doppler factor, the 
magnetic field strength, and the angular size of the optically thick core were derived using published X-ray data. Analysis of 
5 GHz VLBI data reveals the existence of further jet components within the central 5 mas. Additionally the data published so 
far on the GHz-spectrum were discussed at all angular resolutions. Jl 101+7225 turns out to be a standard quasar for studying 
different aspects of radio jet kinematics out to kpc-scales. 

Key words. Galaxies: active - Galaxies: jets - Galaxies: kinematics and dynamics - quasars: individual: Jl 101+7225 - Radiation 
mechanisms: non-thermal - X-rays: galaxies 



1. Introduction 

The radio-loud quasar J 11 1 +7225 is a source of the Very Long 
Baseline Array (VLBA) CaUbrator Survey VCSl for phase- 
referencing observations (Beasley et al. 2002, in the follow- 
ing B02). We observed this source within larger experiments, 
that were focused on imaging faint Seyfert galaxies. It turned 
out that J 1 101+7225 shows a complex source structure, which 
makes it physically interesting but, on the other hand, also less 
suitable as a calibrator. Since so far very little information is 
available for individual sources of the VCSl, we present our 
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property 


value 


ref. 


lAU " / other name: 


J1101+7225/[HB89] 1058+726 




R.A. (J2000) 


llh 01m (48.8054 + 0.0001)s 


B02 


Decl. (J2000) 


+72°25'(37. 1183 +0.6- 10-')" 


B02 


redshift 


1.46 


(b) 


abs. magnitude 


-27.3 mag« lO*" * erg s"' 


(c) 


1 .4 GHz radio cont 


(1451 ±30)mJy 


(d) 



Table 1. General properties of J 1101+7225. ": International 
Astronomical Union; (b): Jackson & Browne (1991); (c): Veron-Cetty 
& Veron (2001); (d): White & Becker (1992) 



results in this article. Analysis of the (relativistic) kinematics 
of non-thermal radio sources gives important physical insight 
into the inmost regions of an active galactic nucleus (AGN). 
The ejection of radio jet components and the jet kinematics are 
most probably related to the process of accretion onto the nu- 
cleus itself. Quantitative estimates of physical properties, such 
as magnetic fields and source sizes, can be derived, too. Our 
analysis is based on comparison of our 1.66 GHz Very Long 
Baseline Interferometry (VLBI) map obtained in 2002 with the 
available maps of past observations (Thakkar et al. 1995 (T95) 
at 1.66 GHz and B02 at 2.3 GHz). Further data from a multi- 
epoch study at 5 GHz of the central 5 mas of J 1101+7225 by 
Britzen (2002) and Britzen et al. (in prep.) were used to deter- 
mine the source structure of the central region, unresolved at 
lower observing frequencies. 

In Table 1 the general properties of J 1101+7225 are sum- 
marized to introduce the quasar B02 estimated the position of 
Jl 101+7225 with (sub)mas-accuracy (Table 1), which is nec- 
essary for high-sensitivity phase-referencing experiments. In 
Sect. 2 details are given concerning the VLBI observations 
that are the basis for this article. Flux densities, the resulting 
maps, and a kinematical analysis of the mapped structure are 
presented in Sect. 3. Finally the results are discussed and sum- 
marized in Sect. 4. 
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date of observation 
13 Feb 2002 



obs. frequency 
1630.49MHz 



bit rate & type of samp. 
256 Mbps & 2 bit 



No. & bandwidth of IF 
4 DSB-IFs & 16 MHz 



polarization 
LCP 



system and correlator: 
MKIV at MPIfR, Bonn 



Table 2. Details of the observing schedule of the EVN observation 
with eight antennae of Jl 101+7225. 



2. Observations 

The quasar Jl 101 +7225 was observed as a calibrator source 
within a larger European VLBI Network (EVN) experiment. 
The observations were conducted during the period 13-14 Feb 
2002, 21h-09h UT, while the scans of Jl 101 +7225 typically 
lasted for three minutes at a total observing time of 1.6 h. 
We observed at 18 cm with the 100m antenna of the Max- 
Planck Institut fiir Radioastronomie (MPIfR) at EfFelsberg, 
Germany, the 76m antenna at Jodrell Bank, UK, the 32m an- 
tenna at Medicina, Italy, the 25m antenna at Onsala, Sweden, 
the 25m antennae at Shanghai and Urumqi, China, the 32m 
antenna at Torun, Poland, and the 14x25m antenna array at 
Westerbork, the Netherlands. Details concerning the observ- 
ing mode are given in Table 2. After the observations, the 
data were correlated at the VLBI correlator of the MPIfR in 
Bonn, Germany, and imported into the Astronomical Image 
Processing System (AIPS) via MK4IN (Alef & Graham 2002). 
The data were fringe-fitted and calibrated in a standard manner 
with the AIPS package and imaged using the DIFMAP VLBI 
package (Pearson et al. 1994). 

Further we present 5 GHz VLBI maps of Jl 101+7225, 
obtained in VLBA and global VLBI observations as part 
of a multi-epoch VLBI study (the Caltech- Jodrell Bank flat- 
spectrum sample CJF; Taylor et al. 1996; Britzen 2002). These 
observations aim at a statistical investigation of the kinematics 
of a complete sample of AGN (Britzen 2002; Britzen et al. in 
prep.). The sources were observed in 5.5 minute snapshot ob- 
servations to determine the position and motion of jet compo- 
nents in the central 5 mas. The data were recorded over 32 MHz 
total bandwidth broken up into 4 baseband channels with 1-bit 
sampling. The recorded data were correlated in Socorro, USA. 

In Sect. 3.2 we include the VLBI-maps of Jl 101+7225 
published earlier to analyze the evolution of the structure 
found. T95 observed the quasar with global VLBI at 1.66 GHz 
in the CJl survey which later on became a part of the CJF sam- 
ple. The other map at 2.3 GHz was reduced from VLBA data 
(B02) via automatic imaging using the Caltech DIFMAP pack- 
age. 

3. Radio properties of the quasar J1101+7225 

3.1. Flux densities 

The flux density of Jl 101+7225 has been measured occasion- 
aUy during the last two decades. In Table 3 we summarize the 
measurements at the various frequencies. Based on these data 
we calculated the spectrum, shown in Fig. 1. The single-dish 
observations show a curved radio spectrum with a steep spec- 



obs. type 


freq. [GHz] 


flux dens. [mJy] 






0.038 


n 41 . 10"^ + SOO"! 




single-dish 


0.178 


(3900 ± 600) 


(b) 


single-dish 


1.4 


(1451 ± 30) 


(c) 


single-dish 


2.7 


(1070 ± 35) 


(d) 


sinclc-dish 


5 


(858 ± 76) 


(c) 


singlc-dish 


12 


(820 ± 100) 


(D 


local interf. 


1.4 


(748 ± 40)/beam 


(f) rcs~l 5" 


local interf. 


8.4 


(349 ± 20)/beam 


(h), res-0.2" 


VLBI 


1.66 


(396+20) 


EVN Feb 2002 


VLBI 


2.3 


(520 + 25")/beam 


Rn2 res-^S mas 


VLBI 


5 


(139 ±7) 


- 1991.4 - 


VLBI 


5 


(282 ± 14) 


- 1993.4 - 


VLBI 


5 


(337 ± 17) 


- 1996.6 - 


VLBI 


5 


(439 ± 22) 


- 1999.9 - 


VLBI 


8.4 


(383 ± 20)/beam 


B02,res~l mas 



Table 3. Measured radio flux densities or peak brightnesses of 
Jl 101+7225, obtained with single-dish telescopes, local interferom- 
eters, and VLBI. The VLBI flux densities are our fitted core flux den- 
sities as presented in Tables 5 & 7 

References and comments: (a) Rees (1990), catalogue revised by 
Hales et. al (1995); (b) Gower et al. (1967); (c) White & Becker(1992); 
(d) Kuhr et al. (1981); (e) Gregory et al. (1996); (f) Terasranta et 
al. (2001); (g) Xu et al.(1995) gave this peak brightness of a map, 
which was convolved to a circular beam of 1.5" FWHM; (h) Patnaik 
et al. (1992) found this peak brightness to be 80% of the total flux 
density of 436 mJy; 

The values are plotted over frequency in Fig. 1 . For the interferometric 
observations 5% flux density errors are given. 




flux densities of Jl 101 +7225 



: single-dish 
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Fig. 1. The radio spectrum of Jl 101+7225 at different angular resolu- 
tions, corresponding to Table 3. All VLBI measurements are labelled 

with the observing date. While the flux densities are plotted with er- 
ror bars, the stars indicate peak brightnesses. Because the peak bright- 
nesses may underestimate the flux densities, the corresponding spectra 
are dotted / dashed (see Sect. 3.1). 



considerable flattening towards higher frequencies (a'-^^f*^ 



-0.4; a\ 



5G;22G 
single 



0). The spectral trend of the Very Large Array 



trum at low MHz frequencies {a\ 



,38M;178M 
single 



-0.8; 5 ~ v") and 



(VLA) measurements follows the single dish measurements. 

The slightly flatter shape {a\^* '^'^ 0.3) plotted in Fig. 1 

may overestimate the real spectrum, because only the peak 
brightness is published at 1.4 GHz. 
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SsGHzLmJy] 


date of observation 


ref. 


ann.: 


778±4 


between Feb 1977 & Mar 1978 


(a) 


(1) 


788±31 


Nov 1986 


(b) 


(2) 


953±73 


Oct 1987 


(c) 


(2) 


623±56 


1992.5 = Jul 1992 


(d) 


(3) 



Table 4. The published single-dish flux densities at 5 GHz to inves- 
tigate the source variability. (1) observed at 4.9 GHz and corrected to 
5 GHz via a spectral index Q-j^q™'; (2) at 4.85 GHz; (3) at 4.75 GHz; 
(a): Kiihr et al. (1981); (b): Gregory et al. 1996); (c): Gregory & 
Condon (1991); (d): Reich et al.(2000) 



In contrast with the EVN we measured a core flux den- 
sity of (396 ± 20) mJy at 1.66 GHz. In combination with the 
core flux densities obtained at 5 GHz and the lower flux den- 
sity limits given by the peak brightnesses of B02, this confirms 
a flat spectral shape (within the scatter due to different obser- 
vations and epochs) at GHz-frequencies, which is shown as a 
dotted line in Fig. 1. These findings indicate that the contin- 
uous flattening of the single-dish spectrum is evoked by su- 
perposition of spectrally steep components, still unresolved 
with local interferometry but fully resolved and/or separable 
at VLBI resolution, and by the flat VLBI core, which is still 
nearly unresolved at the observing frequencies of 1.66 GHz 
and 2.3 GHz. Konigl (1981) calculated such flat GHz-spectra 
of compact synchrotron sources, involving optical thickness for 
the synchrotron radiation due to synchrotron self-absorption. 
Thus most of the central radio emission of Jl 101+7225 at these 
frequencies is emitted by a compact, unresolved core. 

A significant fraction (up to 50%) of the single-dish flux 
density of the whole galaxy is radiated by the inmost region, 
still unresolved by interferometric observations. Relativistic 
beaming models explain the extraordinary luminosity of radio- 
loud cores with unresolved radio-jet components approach- 
ing nearly along the line of sight at relativistic velocities (e.g. 
Blandford & Konigl 1979, and see Sect. 3.4). 

One consequence of these models is a straightforward ex- 
planation of the often observed variability of the core flux den- 
sity by small variations in the jet orientation with respect to 
the observer. This variability is still observable with single-dish 
observations in case of core-dominated GHz emission. Table 4 
presents single-dish radio flux densities, measured at 5 GHz 
at different epochs. The values show a flux density variation 
up to 20% between 1986 and 1987. This variability has to be 
taken into account and leads to significant changes of short- 
range' spectral indices in the GHz domain from different ob- 
serving epochs. The underlying variability of the VLBI core 
flux density was observed at 5 GHz (Table 3) and shows the 
same minimum at around 1992. Therefore only the later mea- 
surements were taken in Fig. 1 to approximate the spectral in- 
dex of the VLBI core. The observed flux density variability of 
the core component strongly supports the concept of beamed 
radio emission, which can be confirmed in Sect. 3.4. 



' Due to the power-law relation between flux density and frequency, 
the influence of the flux density variabiUty on the spectral index de- 
creases with increasing frequency interval. 
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Right Ascension (mas) 



Fig. 2. Cleaned map of the EVN observation of Jl 101+7225 at 
1.66 GHz. The map shows a peak brightness of 407 mjy/beam. The 

contour levels are -0.15, 0.15, 0.3, 0.6, 1.2, 2.4, 4.8, 9.6, 19.2, 38.4, 
76.8% of the peak brightness, and the beam size is 6.22x5.28 mas 
at position angle (PA) of 4°. In the upper right comer the Gaussian 
model fit representation of the VLBI image is given (with the same 
contour levels, but now with respect to a 399 mJy/beam peak bright- 
ness), obtained by fitting Gaussian components to the visibiUty data 
(beam: 5.5x4.3 mas at PA 2°). By this technique the extended struc- 
ture at ~ 50 mas distance to the core could be clearly resolved into 
two components. 

3.2. High apparent superluminal motion outside tlie 
central 10 mas 

In Fig. 2 we present the map made from the EVN observation 
in Feb. 2002. Clearly two extended regions are visible beside 
the core with more than 4.8% and 1 .2% of the core peak bright- 
ness. They are labeled with A and B, respectively. More infor- 
mation about the underlying source structure can be retrieved 
by fitting circular and elhptical Gaussians to the visibility data 
and optimizing the fit to the amphtudes and closure phases (cf. 
Pearson 1995). The structure in 200 mas distance to the core in 
N-NE direction could be confirmed neither by applying differ- 
ent taper to the uv data nor by a respective model component. 
Our best model of the source structure is given in Table 5 and 
the corresponding map is inserted in Fig. 2. The measurement 
errors on the model parameters flux density, relative position, 
and source size are given in the respective tables. Their esti- 
mation is based on the comparison of different model fits, as 
described by Gomez & Marscher (2000). A more theoretical 



4 J.-U. Pott et al.: Relativistic jet motion in the core of the radio-loud quasar Jl 101+7225 



epoch & observ. freq. [GHz]: 


25 Sep 1991 at 1.66 [T95] 


13 Feb 2002 at 1.66 


mean app. transv. velocity 


19 Apr 1995 at 2.3 


components: 


Distances in Fig 


. 3 [mas] 


/i[masyr-i] j8app;*=o.7i 


Dist. in Fig. 3 [mas] 


Core-A 


(19.6 ± 0.5) 


(23.3 ±0.5) 


(0.36 ± 0.07) (22.5 ± 4) 


(21.5 ±3) 


Core-B 


(53.7 ± 1) 


(52.1 + 1) 


0(«) 


(52 ± 5) 



Table 6. The angular distances between the core and the fitted Gaussian jet components (or their flux density weighted mean; cf. text). The 
errors reflect uncertainties in the fitting process and increase with decreasing component flux density, while derivation of the apparent velocities 
is described in the text. The position angles of the components did not change significantly during the observing epochs and were omitted. In 
the last column the read out values from a map by B02 are given for comparison, but cannot be included in the calculations due to lack of a 
Gaussian model fit. 



'"': No significant propagation of the B component was found with respect to the errors. 

*^: For the calculations we used Hq = 71*3 km s"' Mpc"' and a deceleration parameter = 0.1. 





Flux dens. 


location 




size 




(mJy) 


r (mas) 


& (deg) 


a (mas) 


Axial ratio 


Core 


396 








2.88 


0.41@9.4° 


C 


149 


3.3 


14.1 


1.83 


0.39@-89.4° 


A 


54 


23.3 


21.0 


3.62 


1 


Bl 


16 


47.4 


17.7 


6.22 


1 


B2 


7 


62.8 


18.5 


4.21 


1 



Table 5. The best fit Gaussians as mentioned in the text, a labels the 

major axis of elliptical Gaussians. A, Bl, and B2 were only fitted by 
circular Gaussians, because the locations of components of lower flux 
density can be fitted more reliably with less degrees of freedom. The 
flux density errors are ~ 5%, the uncertainties of the location, and size 
of the Gaussians range between 0.2-0.5 of the beamsize depending on 
the respective flux density. 



approacli (Fomalont 1989) based on tlie dynamic range of the 
map gives comparable magnitudes. 

This analysis shows that the extended B component from 
the cleaned map is a blend of at least two components (Bl, 
B2). Furthermore within the central 5 mas a significant fraction 
of the flux density is radiated by an extended source beside the 
dominating unresolved core at a distance of a few mas to the 
core (see Sect. 3.3). This region close to the limit of resolution 
of the 1.66 GHz observations is studied in more detail in the 
next section (3.3) at higher frequencies. 

T95 detected both A and B, too. A detailed comparison of 
our data with the Gaussian component model of T95 reveals 
a high separation velocity of the A-component at an unusually 
large deprojected distance to the core. In contrast to Table 5, 
they could fit three elliptical Gaussian components to the region 
around A and one to the B region. Because the angular exten- 
sions of both regions (A and B) are close to the limit of resolu- 
tion, we prefer the blend interpretation of this situation. Our A 
component does not correspond to one of the three fitted com- 
ponents of T95 but does correspond to a flux density-weighted 
mean of all three components; and vice versa, the mean of our 
Bl and B2 components corresponds to the one component of 
that region fitted by T95 to their visibihty data. Although the 
full width at half maximum (FWHM) of the circular Gaussians 
Bl and B2 (cf . Table 5) are not overlapping, the modehng pro- 
cess showed that due to their low flux densities, the sizes of 
the B components are insecure up to a factor of two. It is also 
possible that a third component exists, slightly too weak to be 



T t -X \ , , , I I I I t I \ r 



T{l,l.660Hz) Tg^JGHz) TP,l.€6GIfe) 




J V 1 1 1 1 I 1 1 

observing epoch 



Fig. 3. The cleaned maps are aligned vertically (rotation about -20°), 
which is justified by the fact, that no significant change of the position 
angle of the jet axis was observed. The apparent linear velocities of the 
fitted Gaussian components (cf. Table 6) are given. The first map by 
T95 at T(l,1.66 GHz) shows a peak brightness of 476 mJy/beam, the 
contour levels are -2,-1, 1, 2, 4, 8, 16, 32, 63, 127, 253 mJy/beam, and 
the beam size is 4.2x3.2 mas at PA 90°. The second map by B02 at 
T(2,2.3 GHz) shows a peak brightness of 520 mJy/beam, the contour 
levels are -3, 3, 6, 12, 24, 48, 96, 192, 384 mJy/beam, and the beam 
size is 3.3x6.7 mas at PA 20°. The beam PAs given here refer to the 
unrotated coordinate system. 



fitted. Thus the used flux density-weighted mean appears to be 
the most adequate presentation. 

We determined, therefore, the apparent motion of the ra- 
dio jet structures A and B with respect to the core (Table 6), 
where in case of several fit components the labels refer to the 
flux density-weighted mean of these. Angular velocities n are 
presented in Table 6 on the right side. They are confirmed by 
the snapshot VLBA map at 2.3 GHz of B02 (last column in 
Table 6). We did not directly include these data in the velocity 
calculation. The introduced uncertainties due to the difl^erent 
observing frequency and the lack of a Gaussian modelfit would 
annihilate the increased calculation accuracy. To demonstrate 
the apparent angular motion, we show the three VLBl maps in 
one figure in Fig. 3 on a common angular scale and rotated by 
about -20 °. 
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5 



epoch 




Flux density 


location 


size 






(mJyj 


r (nuis) 




diam (mas) 


1991.4 


Core 


139 








0.57 




CI 


132 


1 .o 1 


n 9 
u.z 


1.01 




C2 


133 


^ ^8 
o.oo 


L L.J 


0.85 


1993.4 


Core 


282 








0.51 




CI 


113 


7 4.7 


A 


0.57 




C2 


61 


3.76 


10.7 


0.20 


1996.6 


Core 


337 








0.40 




CI 


97 


2.32 


6.0 


1.13 




C2 


80 


3.9 


12.0 


0.94 


1999.9 


Core 


439 








0.41 




CI 


78 


2.07 


7.8 


1.00 




C2 


84 


3.83 


11.1 


1.18 



1058+726 at 4.992 GHz 1991 Jun 09 



Table 7. The best fit circular Gaussians of the 5 GHz data over four 
epochs. The flux density errors are ~ 5%, and the uncertainties of the 
location and size of the Gaussians range between 0.2-0.5 of the beam 
size depending on the respective flux density. 



The dimensionless linear equivalents ySapp are calculated via 



app 



1 + ^/T+2qoz + z 



Ho(l + Z) [ 1 + ^Jl+2qoz + qoz. 



(cf. Pearson & Zensus 1987). The measured apparent superlu- 
minal transverse motion can be transformed via speciid relativ- 
ity^ into a minimal intrinsic velocity, expressed as: 



jSi.„ - 0.9990 



+0.0003 
-0.0005 



nt M X 2 54° ^"-^^ 

^min ~ ^■-'^ -0.38 



(1) 



The 22 mas distance to the core at = 2.5° can be deprojected 
to a large Unear distance of about 4 kpc (atz = 1 .46, = 0. 1). 
It turns out that the uncertainties in the estimation of jSapp and 
the Hubble parameter are affecting by less than 0. 1%. The 
more distant B-component was found to be stationary within 
the errors at a mean angular distance of 52.9 mas. 

3.3. Resolved structure in the central 10 mas 

The higher-resolution maps of the 5 GHz snapshots of Britzen 
et al. (in prep.) are presented in Fig. 4. They explicitly show 
that the Core region of the 1.66 GHz map is not totally com- 
pact on the 1 mas scale of the 5 GHz observations. Extended 
Gaussian components could be fitted to the data beside the 
dominant^ central source. This was already suggested by the 
C-component, fitted to the 1.66 GHz data (Table 5), and it con- 
firms the power of the model fits in analysing the data. 

In Table 7 the fitted Gaussian components are shown. A lin- 
ear fit of the core-distance with respect to the observing epochs 
(Fig. 5) demonstrates the apparently superluminal motion of 
both components and shows mean separation velocities of 



^„j,p^i2.5±2.5) and 



nC2 



app 



(3.7 ± 2.5) 



(2) 



If & is the angle enclosed by the line of sight and the direction of 



motion and if fS is the intrinsic velocity, one finds ySa, 



P sin ,•) 
l-p cos# 



The 



minimal intrinsic velocity jSmin = -^(/Sapp)/(1 filpp) implies an angle 

l?nmi fulfilling: C0t(l?mm) = fiapp- 

^ e.g. at the last epoch the extended components show flux densities 
of 20% of the central flux density 




1058+726 at 4.992 GHz 1993 Jun 10 




Core 



Right Ascension (mas) 
1058+725 at 4.991 GHz 1996 Aug 20 





Right Ascension (mas) 



Riglit Ascension (mas) 



Fig. 4. 5 GHz maps of the central 5 mas. The details are given in the 
form: peals brightness; contour levels in %; beam size. First epoch: 
0.097 mJy/beam; levels are -8, 8, 12, 18, 27, 40.5, 60.8 %; beam: 
1.08x0.79 mas at PA -70°. Second epoch: 0.209 mJyA)eam; levels are 
6, 9, 13.5, 20.3, 30.4, 45.6, 68.3 %; beam: 0.96x0.90 mas at PA -90°. 
Third epoch: 0.314 mJy/beam; levels are 1.6, 2.56, 4.1, 6.55, 10.5, 
16.8, 26.8, 42.9, 68.7, %; beam: 1.95x1.64 mas at PA -67°. Fourth 
epoch: 0.417 mJy/beam; levels are -0.9, 0.9, 1.44, 2.3, 3.69, 5.9, 9.44, 
15.1, 24.2, 38.7, 61.8 %; beam: 2.11x1.64 mas at PA -65°. Crosses 
indicate the positions of circular Gaussian modelfit components. The 
detailed model parameters are presented in Table 7. 



The given standard deviations of the fitted velocities are rel- 
atively large due to the few available data points. Because a sig- 
nificant intrinsic acceleration of the jet material is very unlikely 
from the 4 mas regime of the inner C-components toward the 
20 mas regime far away from the central engine, we adopt the 
yg^.^ ss 0.999 from Eq. 1 as the intrinsic separation velocity also 
for the inner jet components. From the observed apparent sepa- 
ration velocities of CI and C2 one uses the equations of special 
relativity ^ to calculate two possible jet orientation angles ■fisj to 
the line-of-sight for each component (with jS^^ = 0.999+o o°°3): 



j8£! * (2.5 + 2.5) : e [0°, 0.45°], e [22°, 180°] 



app 
qC2 



f^app ~ (3.7 + 2.5) : £ [0.06°, 0.56°], e [18°, 80°] (3) 



J.-U. Pott et al.: Relativistic jet motion in the core of the radio-loud quasar Jl 101+7225 



C2: 0.06 mos/yr 





CI: 0.04 mos/yr 



radio 


X-ray 


results 


5m@Vm; V2 


^(lkeV;5keV) „) 


Doppler factor 


300±15mJy@5 GHz; 2 THz 


-0.5 + 0.2 


S = 0.5:j;^ 


^irc 




magnetic field fio 


(0.9 ± 0.45) mas 


0.1//Jy 


(0.15t?2)G 



Table 8. The radio values characterize the compact core at 5 GHz as 
discussed in Sect. 3.4. The turn-over properties are mean values of 
the model fits in Table 7. The cut-off frequency V2 of the synchrotron 
spectrum is extrapolated from the optical-radio spectrum. The right 
column gives the calculated values. From Fiore et al. (2001; 5 keV) 
and Brinkmann et al. (1997; 1 keV). 



Fig. 5. Evolving core-distance (in mas) of the best fit Gaussians, 
taken from Table 7 over the years. A least-squares fit of a linear 

time-dependence is shown with respect to the errorbars. The stan- 
dard deviation of the resulting apparent separation velocities is cr « 
0.04mas yr"^ 



The uncertainty of ySfpp^ dominates the errors of Eq. 3. For both 
components the calculated orientations differ significantly from 
<in - (2-54!°J^r ofEq.l. 

A closer look at the modelfits in Table 7 reveals some 
trends over the observing epochs. The separation motion seems 
to decrease at later epochs in correlation with a decrease in 
the respective component brightnesses and an increasing posi- 
tion angle of CI from 0.2° to 7.8° over the four epochs. With 
the given, relatively poor time-resolution of four observations 
within eight years it is not possible to fit these trends reliably. 
Nevertheless in combination with the different calculated an- 
gles in Eq. 1&3, they can be interpreted as indications of a spa- 
tially curved jet structure. Different line-of-sight orientations 
account for different brightnesses, position angles, and appar- 
ent separation velocities. 

The core flux density has strongly increased over the four 
epochs, which suggests that over the later epochs either a new 
jet component has emerged from the compact core region or 
the unresolved jet has changed its orientation with respect to 
the line of sight. During the fitting process we found indications 
that indeed a third component in sub-mas distance to the core is 
hidden; but probably due to its proximity to the core, this new 
component could not be fitted separately without doubts. 

On the other hand, the increased brightness of the core 
at the later epochs may affect uncertainties in the model fit- 
ting process of the close C components. The values of one 
fifth of the beam size (cf. caption of Table 7 and discussion 
in Sect. 3.2), which are used for the linear fit in Fig. 5, may 
still shghtly underestimate the real errors due to increased 
overblending of the C components by the close core radiation. 
Because this effect is not quantifiable with our data, we kept 
using these errors. 

The idea of a jet curvature over the observed core-distances 
does not rule out a ballistic situation at the origin of the jet, as 
was described recently by Stirling et al. (2003). As origin of 
the radio jet in BL Lacertae they found a 'precessing nozzle' 



which ejects the single components along straight trajectories. 
But outside 2 mas (corresponding to 0.4 mas at the redshift of 
JllOl-l-7225) these trajectories became curved as well. 

Summarizing the motion of all detected jet components in- 
cluding A and B, the situation resembles the different measured 
apparent component speeds of the radio jet of the S5 quasar 
0836-H710 which extends over more than 150 mas at z ~ 2.17 
(Hummel et al. 1992) and where no systematic correlation be- 
tween the component speed and its distance to the core seems 
to be present (Otterbein et al. 1998 and references therein). 

3.4. Indirect estimation of relativistic bulk motion in the 
core 

Our VLBI-data of Jl 101-1-7225 allow, in combination with pub- 
lished data at other wavelengths, a straightforward analysis of 
the observed VLBl core flux density by applying a relativistic 
beaming model. The crucial idea of such models is that due to 
relativistic bulk motion of the sources, the radiation is ampli- 
fied or attenuated significantly in the rest-frame of the observer, 
depending on the direction of motion with respect to the Une of 
sight. Of course the bulk motion cannot be observed directly 
due to the lack of spatial resolution. But we share the common 
assumption that the X-ray emission of the core consists mainly 
of synchrotron photons scattered to shorter wavelengths via the 
inverse Compton effect. Marscher (1983) quantitatively derived 
the connection between radio and X-ray flux densities, their 
spectral indices, and the Doppler factor"* 6 for a homogeneous 
source of spherical shape. This implies that 6 can be estimated 
from the other measured values. 

In Table 8 the ingoing values of the calculations are given. 
The notation of the radio properties follows Marscher (1983). 
Typically the unresolved core consists of several components 
each of which dominates the unresolved spectrum at a different 
frequency (e.g. Marscher 1988). Therefore we can assume that 
the measured properties at 5 GHz describe the turn-over point 
(index m) of one underlying source component and that the 
formulae can be apphed for that component. The X-ray mea- 
surements were taken independently in the 1990s. Therefore a 
mean of the different 5 GHz measurements is the most appro- 
priate way to combine both datasets. 



** Wit h the te rms of footnote litis 6 = (y(l -yS cos !?))"' 

r= 1/Vi^ 



with 



J.-U. Pott et al: Relativistic jet motion in the core of the radio-loud quasar Jl 101+7225 



7 



The mean measured Gaussian source size is transformed 
into 0'"'^ assuming an underlying circular source shape (see 
Marscher 1983 and Table 8). The angular source size should be 
understood as an order of magnitude estimate^. Further we can 
use only the measured peak flux density 5^ in contrast to the 
extrapolated one as foreseen by Marscher (1983). This results 
in a slight underestimation of 6. The calculations lead to 5 = 
(0.5:^Jj) and a magnetic field strength of Bo=iOA5t.°Qj^) G of 
the unresolved synchrotron-self absorbed source (cf. Table 8). 

However further radio components, peaking at higher fre- 
quencies, and other processes than inverse Compton scattering 
may contribute to the observed X-ray flux densities. This re- 
sults in a probable underestimation of 6 and Bq. In fact a flat 
to inverted spectral index as indicator for several components 
is suggested by the VCSl map at 8.4 GHz (B02). Although the 
lack of a detailed model fit inhibits reproducible calculations, 
visual inspection of the map reveals a dominating core, prob- 
ably unresolved. Under the latter assumption the given peak 
brightness resembles the flux density of the component in mid 
1995, which is slightly larger than the respective values of the 
5 GHz experiment at that time (cf. Fig. 1). 

With the angular source size and the Doppler factor, the 
intrinsic brightness temperature can be calculated as 

Ti„,r = T,pp/d^(2.Slli-10''K) . (4) 

These values fit well in the range of brightness temperature in 
which the simultaneously observed synchrotron emission and 
the inverse Compton-scattering are both effective enough to 
fulfill the inverse Compton scenario adopted here (2 • lO" K < 
Tic < 10^^ K; e.g. Kraus 1986, Bloom & Marscher 1991). 
Without applying a Doppler factor the apparent brightness 
temperature at 5 GHz rises from 0.5-3.1-10^' K over the four 
epochs. Both this strong rise and the fact that the earlier esti- 
mates are below the lower inverse Compton linoit can be inter- 
preted as indicators of relativistic beaming. 

Furthermore we calculated the equipartition Doppler factor 
5equ following Readhead (1994) and Guijosa & Daly (1996), 
which implies that the radiating particles have the same en- 
ergy as the penetrating magnetic field. The results (assuming 
as above = 0.71) 

Se,u = (0.5t_li) (5) 

are very similar to the inverse-Compton Doppler factor 
(Table 8). This supports Readhead's conclusion that many pow- 
erful, non-thermal extragalactic radio sources are close to en- 
ergy equipartition and their finding that the brightness temper- 
atures of the respective sources range far below the maximum 
brightness temperature of lO'^ K (see above). 

As in the previous section, we now adopt that the intrinsic 
velocities along the jet are equal to the highest minimal intrin- 
sic velocity, as given by the apparent separation speeds. Then 

' If much more (spectral and multi-epoch) data are available as e.g. 
in the case of the quasar 3C 345, inhomogeneities in the particle num- 
ber density can also be estimated. Lobanov & Zensus (1999) used 
more detailed analyzis to explain the observed flux density variations 
quantitatively with time of 3C 345. 



the angle between the direction of motion and the line of sight 
can be estimated from the Doppler factor: 

Score = iOXliyXin ^ (0.999!°;°°°^) ^ ^^^^^ « (24:J^)° (6) 

Thus ■d'coie includes only the large angle solutions of Eq. 3 and 
supports the idea of a straight inner jet nozzle as described in 
Stirling et al. (2003). 

4. Discussion of the results 

For the first time the kinematics of the non-thermjil radio jet 
components of Jl 101+7225 were estimated. We calculated the 
apparent motions of the two extended components, observable 
at low GHz-frequencies from VLBI-measurements. We found 
an apparent superluminal separation speed of /3app(Core -A) - 
(22.5 ± 4) for the A-component over the past decade at an ex- 
ceptionaUy large deprojected distance to the core (22 mas ~ 

4 kpc). The even more distant B-component was found to be 
stationary within the errors at a core-distance of ~ 53 mas. 

Furthermore a Doppler factor S = (0.5;|;g 3) was estimated 
for the compact optically thick (at the observing frequencies) 
synchrotron radiation of the core of J 1 101+7225. The strength 
of the magnetic field, which is necessary for the synchrotron 
process, was estimated to be Bo = (0.15^^^^) G. The esti- 
mated size of the optically thick nuclear component is close 
to the resolution limit of an interferometric observation at cm- 
wavelengths, which suggests, in combination with the esti- 
mated Doppler factor and intrinsic jet velocities, that additional 
VLBI-components could appear outside the core in the future. 

More observations of these nuclear components could re- 
veal, if a deceleration of the A-component appears while A 
is approaching the B-component. This would be similar to 
file well known case of 4C 39.25 (Alberdi et al. 1993a). 
Deceleration at these distances from the core can confirm an in- 
teraction of the jet component with the circumnuclear matter of 
the host galaxy as described by Taylor et al. (1989) for Seyfert 
nuclei. Such an interaction is strongly supported by compar- 
ing the VLBI maps presented here with the jet geometry as ob- 
served with the VLA (Xu et al. 1995). In their 1 .4 GHz map two 
corresponding radio jets appear, extending along ~ 5" from the 
core. While the southwestern jet extends in the opposite direc- 
tion of the VLBI jet, the northeastern VLA-jet does not coin- 
cide with the position angle of the VLBI structure. Instead it 
appears to be bent towards a northwestemly direction, perhaps 
induced by ram pressure of the surrounding material. 

Apparent superluminal velocities are explained by motion 
towards the observer at relativistic velocity. The additionally 
presented 5 GHz maps and Gaussian modelfits of the central 

5 mas show further jet components with significantly smaller 
apparent velocities. Thus the separation velocities of the differ- 
ent VLBI radio jet components of Jl 101+7225 show no sim- 
ple correlation with their distance to the core. But this does 
not rule out a constant intrinsic velocity, because already smaU 
variations of the angle between jet and the line of sight can in- 
troduce variations of the apparent speeds of the observed order 
of magnitude. In the case of such a constant intrinsic velocity 
along the whole jet, the very high velocity, derived from the 
A-component, has to be chosen. Such high intrinsic velocities 
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explain both the large extension of the radio jet over several 
arcsec as estimated with the VLA and the high luminosity of 
the core. 

We believe that our findings may indicate a helical bending 
of the jet, where the fast components are moving in a section 
that is curved towards the observer. This is supported by the 
estimated differing jet orientations with respect to the line of 
sight. Zensus et al. (1995) could explain the acceleration of a 
jet component of 3C 345 by a curved jet of constant intrinsic 
bulk velocity. It has been shown that helical jet patterns result 
from Kelvin-Helmholtz and current-driven jet instabilities in 
relativistic flows (Birkinshaw 1991, Istomin & Pariev 1996). 

Beside this interpretation of the observed motion as an in- 
trinsic bulk motion of the radiating plasma, shock waves may 
travel along the jet (Alberdi et al. 1993b). Diff'erent component 
velocities may also be observed if the slower components are 
trailing in the wake of faster ones. This hydrodynamical expla- 
nation could be adopted successfully to the complex compo- 
nent motion in the radio jet of the radio galaxy 3C 120 (Gomez 
et al. 2001). Hardee et al. (2001) found mechanisms to produce 
differentially moving and stationary features in a jet by analyz- 
ing the relativistic hydrodynamic equations. 

Thus the nuclear region of Jl 101-1-7225 provides the rare 
possibility of observing the total range of jet kinematics includ- 
ing apparent superluminal separation velocities even far out of 
the central parsec-region. The results presented here can give 
observational constraints far from the jet origin for numerical 
jet models. 
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